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Abstract. Starting from the three-dimensional Gross-Pitaevskii equation (3D GPE) 
we derive a ID generalized nonpolynomial Schrodinger equation (ID g-NPSE) which 
describes the dynamics of Bose-Einstein condensates under the action of a generic 
potential in the longitudinal axial direction and of an anisotropic harmonic potential in 
the transverse radial direction. This equation reduces to the familiar ID NPSE [Phys. 
Rev. A 65, 043614 (2002)] in the case of isotropic transverse harmonic confinement. 
In addition we show that if the longitudinal potential models a periodic optical lattice 
the 3D GPE can be mapped into a ID generalized discrete nonpolynomial Schrodinger 
equation (ID g-DNPSE). 
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1. Introduction 

At ultralow temperature dilute Bose-Einstein condensates (BECs) can be accurately 
described by the three-dimensional Gross-Pitaevskii equation (3D GPE) [JJ. Some years 
ago we have found that, starting from the 3D GPE and using a Gaussian variational 
approach [21 [3] , it is possible to derive an effective ID wave equation that describes the 
axial dynamics of a Bose condensate confined in an external potential with cylindrical 
symmetry [HE]. In this derivation the trapping potential is harmonic and isotropic 
in the transverse direction and generic in the axial one. Our equation, that is a time- 
dependent nonpolynomial nonlinear Schrodinger equation (ID NPSE) [HE], has been 
used to model cigar-shaped condensates by many experimental and theoretical groups 
(see for instance [6l[7J[8j[9]). Moreover, by using the ID NPSE we have found analytical 
and numerical solutions for solitons and vortices, which generalize the ones known in the 
literature [10] . Recently, we have obtained a disrete version of NPSE, the ID DNPSE, 
which models Bose-Einstein condensate confined in a combination of a cigar-shaped trap 
and deep optical lattice acting in the axial direction [IT] . 

A relevant limitation of the ID NPSE [U [5] is the fact that the transverse harmonic 
confinement is isotropic. To overcome this problem in this paper we introduce a 
ID generalized nonpolynomial Schrodinger equation (ID g-NPSE) which describes the 
dynamics of matter waves under the action of a generic potential in the longitudinal axial 
direction and of an anisotropic harmonic potential in the transverse radial direction. 
This equation reduces to the familiar ID NPSE [H [5] in the case of isotropic transverse 
harmonic confinement. In the second part of the paper we consider a deep periodic 
optical lattice in the axial direction. In this case we find that the 3D GPE can be 
mapped into a ID generalized discrete nonpolynomial Schrodinger equation (ID g- 
DNPSE), which becomes the ID DNPSE obtained in Ref. [IT] if the matter waves 
are confined by a transverse isotropic harmonic trap. 

2. ID generalized NPSE 

We consider a dilute BEC confined in the axial direction by a generic potential V(z) 
and in the transverse plane by the anisotropic harmonic potential 



where A is the anisotropy of the transverse harmonic trap: oj x = Xuj± and u) y = uj± 
are the two harmonic trapping frequencies in the x and y directions. The characteristic 
harmonic length is a± = (h/(mu±)) 1 ^ 2 and so the characteristic length and time units 
are a±_ and uj 1 , and the characteristic energy unit is hu±. The system is well described 
by the 3D GPE, and in scaled units it reads 




(1) 
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where ip{r,t) is the macroscopic wave function of the condensate normalized to unity, 
and g = 2Na s /a±, with N the number of atoms and a s the s-wave scattering length of 
the inter-atomic potential. This equation can be derived from the Lagrangian density 



- l -\^\ 2 -\{\ 2 x 2 + y 2 



■V{z 



.(3) 



To reduce the 3D problem to a ID one, we apply a variational approach. In the 
case of an anisotropic cigar-shaped geometry, it is natural to extend the variational 
representation [U [5] of the 3D GPE which led to the NPSE for the axial wave function 
in the isotropic case to the present situation. Thus we adopt the following ansatz for 
the state described by 3D GPE 
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where a(z,t), rj(z,t) and f(z,t), which account for transverse widths and axial wave 
function, are the 3 effective fields to be determined variationally. 

By inserting this ansatz into the Lagrangian density ([3]), performing the integration 
over x and y and neglecting spatial derivatives of transverse widths (adiabatic 
approximation [1]), we derive the effective Lagrangian density 

2 
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This effective Lagrangian gives rise to a system of 3 Euler-Lagrange equations, obtained 
by varying C with respect to /*, a and r\: 
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(6) 
(7) 
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We call this system of three coupled equations, which describe the BEC under 
a transverse anisotropic harmonic confinement, the ID generalized nonpolynomial 
Schrodinger equation (ID g-NPSE). Notice that with g ^ only for A = 1 one finds 
a = rj. In this case of isotropic transverse harmonic confinement (A = 1) we have 



a 



V 



g\f\ : 



and the g-NPSE becomes 



l dt f 



1 d 2 



/, 



(9) 



(10) 



^ + g\f\ : 

which is the ID NPSE derived in Ref. [H E]. Instead, for an anisotropic transverse 
harmonic confinement (A / 1) the transverse widths a and rj depend on g in a not 
trivial way. 



Generalized nonpolynomial Schrodinger equations 



4 



2.1. Weak-coupling regime 

In the weak-coupling regime, i.e. g\f\ 2 <C 1, we can expand in series of powers of g\f\* 
the widths a and rj. At the first order in g\f\ 2 , the g-NPSE gives 

1 1 m2 

V = l-^9\f\ 2 
and we get a ID GPE for the field f{z, i) 



.df 

''dt 
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(12) 



(13) 



In this weak-coupling limit the anisotropy produces a renormalization g\f\ of the 
interaction strength g and the transverse energy is simply (A + l)/2. 

2.2. Strong- coupling regime 

In the strong-coupling regime g\f\ 2 3> 1 the g-NPSE gives 

wi 2 ) 1/4 



° A3/4 

n = a 1/4 (^|/| 2 ) 1/4 

and the differential equation reads 
.df 



dt 



\^+nz)\f\ 2 +l^\f\]f 



(14) 
(15) 

(16) 



that is a ID Schrodinger equation with quadratic nonlinearity. In this regime, after 
setting 

f(z,t) = p 1 (z) 1 / 2 exp(-tfit) , (17) 

where \i is the adimensional chemical potential, and neglecting the spatial derivative 
(Thomas- Fermi approximation), one finds the stationary axial density profile of the 
BEC 

4 
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2.3. Quasi-2D regime 

Under the condition g\f\ 2 t] a, which implies a very strong harmonic confinement 

along the y axis (A < 1), Eq. (|Sj) gives r\ ~ 1 and the system is quasi-2D. In this case 
from Eqs. (jSJ) and (JTj) we get 



9 f 
% dt f 



2dz 2 4 \a 2 
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(19) 



and also 



AV = l + g\f\ 2 a 



(20) 
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We observe that this last equation, based on a 3D— >2D— >1D reduction of the GPE, is 
formally equivalent (with A = 1) to an equation for the width a previously obtained 
within a 3D— >2D reduction of the GPE [T2] . Exact solutions to Eq. ( 1201) are given by 
the Cardano formula, 

^ = T^{±^/ £1 ^ (21) 

3A ^"\V^ J 1 ^ 3A 
where the upper and lower signs correspond, respectively, to g > and g < 0, and 

1/3 




A = (3/2) 1 / 3 f^l/l 4 + ^256 + ^|H/| 8 j • (22) 

Analytical-numerical approach 

In general, the solutions of Eqs. ([7]) and (El) do not have a simple analytical expression. 
Nevertheless from Eqs. ([7]) and ([HD one finds that the transverse widths a and r\ can be 
expressed in parametric form as 

^ = 4l/| 2 ( ^ 1)1/4 , (23) 

V = (Z + l) 1/4 (24) 
where the parameter £ is the solution of the quartic equation 

£ 3 (£+s 2 m 4 ) = Aym 8 (£+i). (25) 

The algebraic equation (125]) can be solved numerically. Nevertheless some analytical 
results can be easily obtained in limiting case. Eqs. (123]) and (124]) show that for a 
repulsive BEC {g > 0) it must be £ > 0, while for an attractive BEC (g < 0) it is 
-1 < f < 0. If A = 1 then £ = g\f\ 2 and a = r, = (1 + ^|/| 2 ) 1/4 . If g\f\ 2 = 1 then 
£ = A 2 / 3 and a = n/X 2 / 3 with 77 = (1 + A 2 / 3 ) 1 / 4 . If g\f\ 2 = -1 then £ = -1 and 
a = rj = 0. 

In Fig. [T]we report the transverse widths a and ?y obtained from Eqs. (123|) . ( jMP 
and the numerical solution of Eq. (125]) with the Newton-Raphson method. We plot 
the widths as a function of the axial strength g\f\ 2 for three values of the transverse 
anisotropy A. 

To test the accuracy of ID g-NPSE we consider the case with negative scattering 
length (g < 0) without axial confinement (V(x) = 0). Under these conditions the system 
admits triaxial bright soliton configurations [13j. In Fig. [2] we plot the axial density 
profile \f(z) \ 2 of the BEC with g — — 1 and three values of the transverse anisotropy A. 
The solid curves represent the numerical results obtained by solving the 3D GPE, Eq. 
(J2]), with a finite-difference Crank-Nicolson algorithm with imaginary time [14J. In this 
case the axial density profile is given by 

p(z)= [ mr)\ 2 dxdy, (26) 
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Figure 1. Transverse widths a and r\ as a function of the axial strength g\f\ 2 - Three 
values of the transverse anisotropy: A = 0.1 (dotted curves); A = 1 (dashed curves); 
A = 10 (solid curves). 
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Figure 2. Axial density profile |/(z)| 2 of the triaxial bright soliton, i.e. attractive 
BEC without axial confinement (V(x) — 0). Solid curves: numerical solutions of the 
3D GPE, Eq. ©. Dashed curves: numerical solutions of the ID g-NPSE, Eq. ©. 
Dot-dashed curves: numerical solutions of the ID GPE, Eq. (|13p . Interaction strength 
g = — 1 and three values of the transverse anisotropy A. 
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with ip(r) ground-state wave function of the 3D GPE. The dashed curves are obtained 
by solving the g-NPSE and the dot-dashed curves by the ID GPE. The figure shows 
that the results of g-NPSE are in good agreement with the 3D GPE ones, while the 
predictions of the ID GPE, Eq. (fTBl) . are not very reliable for the interaction strength 
<? = -!■ 



3. ID generalized DNPSE 



Let us consider again Eq. (j2J) where now the axial potential is given by 

V(z) = V cos{2kz) . (27) 

This potential models the optical lattice produced in experiments with Bose-Einstein 
condensates by using counter-propagating laser beams [T5] . 
To symplify the problem we set 

il>(x,y,z,t) = y £<l> j (x,y,t) W,(z) (28) 

3 

where Wj(z) is the Wannier function maximally localized at the j-th minimum of the 
axial periodic potential. This tight-binding ansatz is reliable in the case of a deep optical 
lattice [T5l [16] . We insert this ansatz into Eq. (j2J), multiply the resulting equation by 
W*(z) and integrate over z variable. In this way we get 



■ 9 A 



( \1 2 , 2 

- (A x +y 



-Ji 



J n+1 



+ 



'n-1 
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k, (29) 



where the parameters e, J and 7 are given by 

1 o 2 



w: z) 



2 0z 2 



+ Vo cos (2kz) 



WJz) dz 



j = -J w: +1 (z) 

1 = 9 J \W n (z)\ 4 dz 



1 d 2 



2dz 2 



+ V cos (2kz) 



W n (z) dz 



(30) 

(31) 
(32) 



In the tight-binding regime the parameter J is positive definite. In addition the 
parameters e, J and 7 are independent on the site index n. 

Eq. fl29l) can be seen as the Euler-Lagrange equation of the following Lagrangian 
density 

" d 1 

2 



dt ' 2 ' M 2' 

+ J0* (0 n+ i + n _i) - 7T7|0 n | 4 } . 

To further simplify the problem we set 



(33) 
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where cr n (t), f] n (t) and f n (t), which account for discrete transverse widths and axial 
wave function, are the effective generalized coordinates to be determined variationally. 

We insert this ansatz into the Lagrangian density ( |33i) and integrate over x and y 
variables. In this way we obtain an effective Lagrangian for the fields f n (t) and a n (t). 
This Lagrangian reads 

~; d l ( l xV 4. 1 



(35) 



To obtain this expression we have supposed that the transverse widths a n (t) of nearest- 
neighbor sites are practically equal (this is the equivalent to the adiabatic approximation 
made for the continuous model). The Euler-Lagrange equation of (1351) with respect to 

fn is 



If 

l dt U 



— J (fn+l + /i 



n-1 



+ 



fn 



7 



' I fn | /n • 



while the Euler-Lagrange equations of (1551) with respect to cr n and rj n give 



AV* = l + 7 |/n 



^ = l+7l/n| 2 ^- 



(36) 

(37) 
(38) 



We call this system of coupled equations, which describe the BEC under a transverse 
anisotropic harmonic confinement and an axial optical lattice, the ID generalized and 
discrete nonpolynomial Schrodinger equation NPSE (ID g-DNPSE) In the case of 
isotropic transverse harmonic confinement (A = 1) we have 



<=Vn = ^ + l\fn\ 2 

and the ID g-DNPSE becomes 
.0 



l dt fn 



tfn — J (fn+l + fn-l) + 



( 3/2) 7 |/n| 

+ g\f n \ 2 



' fn 



(39) 



(40) 



which is substantially the ID DNPSE derived in Ref. [TT] . 

The ID g-DNPSE given by Eqs. ([36]), (J37D and (J3SJ is the discrete version of ID 
g-NPSE given by Eqs. (jHJ), (E) and (jHJ). It is then strightforward to repeat with the ID 
g-DNPSE the analysis of weak-coupling, strong-coupling and quasi-2D regimes we have 
performed with the ID g-NPSE. 



4. Conclusions 



We have obtained new nonpolynomial Schrodinger equations for matter waves under 
anisotropic transverse harmonic confinement. In the case of isotropic transverse 
confinement these continuous and discrete equations reduce to the ones derived in 
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previous papers. The generalized ID continuous nonpolynomial Schrodinger equation 
can be used to study the static and dynamics of Bose condensates with a generic 
longitudinal axial potential. The generalized ID discrete nonpolynomial Schrodinger 
equation can be insead used to investigate the effect of a deep axial periodic potential. 
The comparison with the full 3D Gross-Pitaevskii equations show that these effective 
ID equations are indeed quite accurate and much simpler for numerical and analytical 
computations. 

The author thanks Alberto Cetoli, Boris Malomed, Giovanni Mazzarella, Dmitry 
Pelinovsky and Flavio Toigo for useful comments and suggestions. This work has been 
partially supported by Fondazione CARIPARO. 
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